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Selenium impurities in silicon are deep double donors and their optical and electronic properties 
have been recently investigated due to their application for infrared detection. However, a singly- 
ionised selenium donor (Se + ) possesses an electron spin which makes it a potential candidate as 
a silicon-based spin qubit, with significant potential advantages compared to the more commonly 
studied group V donors. Here we study the electron spin relaxation (Ti) and coherence (T 2 ) times 
of Se+ in isotopically purified 28-silicon, and find them to be up to two orders of magnitude longer 
than shallow group V donors at temperatures above ~ 15 K. We further study the dynamics of 
donor-acceptor recombination between selenium and boron, demonstrating that it is possible to 
control the donor charge state through optical excitation of neutral Se°. 


I. INTRODUCTION 

The electron and nuclear spins of Group V donors in 
silicon, such as phosphorus, have been actively studied 
as potential quantum bits mia. Recent advances in this 
field include the manipulation and read-out of individ¬ 
ual electron and nuclear spins by integrating donors into 
nanoelectronic devices EH5], and the demonstration that 
donor spin coherence times can be as long as 3 seconds 
for the electron spin [6j and up to 3 hours for the nuclear 
spin [7]. In contrast to such work on shallow donors in sil¬ 
icon as qubits, there has been relatively little experimen¬ 
tal attention on the chalcogens, which are deep double 
donors in silicon. Nevertheless, they possess several at¬ 
tributes of potential relevance to donor-based spin qubits 
— we focus here on selenium as an example. 

First, neutral selenium (Se°) possesses two bound elec¬ 
trons which form a singlet ground state, and has a bind¬ 
ing energy of 307 meV [5]. This makes it ideal for spin 
readout methods using spin-to-charge conversion that re¬ 
quire that the two-electron state on the donor is well 
bound B mu- In phosphorus, for example, the two- 
electron state (P _ ) has a very weak binding energy of 
only 2 meV. Second, the singly-ionised state, Se + , is iso- 
electronic with the neutral Group V donors and possesses 
a bound 5 = 1/2 electron spin. Compared to shallow 
donors, Se + has a much larger thermal ionisation energy 
(593 meV) [TT] such that it can retain a bound electron 
even at room temperature. Third, the large energy sepa¬ 
ration between the donor ground state, ^(Aq), and first 
excited valley state, ls(T 2 ) is of order 429 meV (more 
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than an order of magnitude greater than for the shallow 
Group V donors) and could be expected to result in a 
significantly longer electron spin relaxation time. 

Here we present a study of the electron spin proper¬ 
ties of Se + in isotopically purified 28-silicon, including 
timescales for electron spin relaxation (Ti) and decoher¬ 
ence (T 2 ), and their mechanisms. We find that the elec¬ 
tron spin relaxation times are as much as two orders of 
magnitude longer than phosphorus, for a given temper¬ 
ature, and that electron spin decoherence times can be 
reasonably expected to be as long as those measured for 
phosphorus (up to seconds [H]). We go on to investi¬ 
gate donor-acceptor (DA) recombination following above 
band gap illumination (1047 nm), by monitoring the elec¬ 
tron spin echo intensity, using the same method previ¬ 
ously applied to phosphorus donors m■ This demon¬ 
strates an optical mechanism for placing the donors in 
the neutral Se° state, thus removing the hyperfine in¬ 
teraction between the electron and the nucleus. We 
find DA recombination is slow (minutes to hours) given 
the concentrations of selenium (~ 5 x 10 15 cm -3 ) and 
boron (~ 5 x 10 13 cm -3 ) used here m and the rate 
at which charge equilibrium is re-established can be en¬ 
hanced by selectively ionising the Se° via optical illumi¬ 
nation at 4 /im. 


II. MATERIALS AND METHODS 

The 28 Si:Se samples used here have been previously 
measured by IR absorption spectroscopy mi, and consist 
of 28 Si doped with selenium, and partially compensated 
with boron in order to produce a significant concentration 
of Se + , as described by Ludwig [15]. The starting silicon 
material had a composition of 99.991% 28 Si, with 75 ppm 
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FIG. 1: X-band cw-ESR of Se+ in 28 Si at X-band for 
28 Si: nat Se (upper) and 28 Si:' 7 Se (lower). The natural abun¬ 
dance of ,7 Se is 7.5%, with the remaining isotopes (92.5% 
abundance) possessing zero nuclear spin. The ,7 Se-doped 
sample shows predominantly the hyperfine-split lines arising 
from coupling to the I = 1/2 nuclear spin of 77 Se, while addi¬ 
tion features in the spectrum correspond to residual isotopes 
of Se and the presence of SeH pairs. Microwave frequency = 
9.38 GHz; temperature = 23 K. 


29 Si and 15 ppm 30 Si, and a residual boron concentration 
[B] ~ 5 x 10 13 cnr 3 . This was then doped with selenium 
by thermal diffusion 0. We studied two samples, one 
with natural abundance selenium, and the other made 
using 77 Se with an isotopic enrichment of 97.1 %. The 
sample resistivity was measured to be 4.4 flcm which cor¬ 
responds to a donor concentration of ~ 5 x 10 15 cm -3 , 
making the samples n-type. ESR measurements were car¬ 
ried out using a Bruker Elexsys 580 spectrometer at X- 
band (9.7 GHz). 



FIG. 2: Temperature dependence of electron spin relaxation 
T\ in 28 Si:' 7 Se + and Si:P at B ~ 0.34 T. Solid lines are fits 
to the experimental data, made up of contributions (dashed 
lines) from different relaxation mechanisms. Data for Si:P is 
from Ref [18]. 


where u> e = g(3Bo/h and u> n = g n /3 n Bo/h are the electron 
and nuclear Zeeman frequencies and B 0 is the static mag¬ 
netic field applied along the z axis. Neutral Se° has a spin 
singlet ground state and thus gives no contribution to the 
ESR signal. Figure [T] shows continuous wave ESR spec¬ 
tra for 28 Si: Ilat Se and 28 Si: 77 Se. The values we obtain for 
A and g confirm previously reported results m , though 
the transitions have a much smaller linewidth (< 5 /iT) 
due to silicon isotopic purification. The relative inten¬ 
sity of the central line (around g = 2, corresponding to 
Se isotopes with zero nuclear spin) and hyperfine-split 
lines (corresponding to 77 Se) in each sample is consistent 
with their expected isotopic composition. The remaining 
features in the 77 Se + ESR spectrum have been charac¬ 
terised by angular dependent cw-ESR and ENDOR and 
confirmed to be due to SeH pairs, previously measured 
in natural Si and referred as Si-NL60 HU- The presence 
of the hydrogen in this defect centre reduces the ' 7 Se hy- 
perfine coupling by a factor of and gives it a slight 
anisotropy. 


IV. ELECTRON SPIN RELAXATION 


III. CONTINUOUS WAVE ESR 

The ESR proprieties of Se + in silicon have been pre¬ 
viously investigated by Grimmeiss et al. [SI SSI- The 
electronic ^-factor is g = 2.0057. All stable Se isotopes 
x Se + (X=74, 76, 78, 80) have zero nuclear spin, apart 
from 77 Se which has I — 1/2 and an isotropic hyper- 
fine coupling of A = 1.6604 GHz with the donor electron 
spin. The electron and nuclear spin proprieties of Se + are 
described by the spin Hamiltonian (in frequency units): 

( 1 ) 


We studied the electron spin relaxation time (Ti) of 
Se + by pulsed ESR, using an inversion-recovery sequence 
(7r — T — 7T/2 — r — 7r — t— echo, where T is varied). The 
results, shown in Figure [2j show that (1 /Tj) is well fit 
by 1/Tj = CT 9 with C = 1.2 x 10 -8 s _1 K 9 , in the tem¬ 
perature range of 5 — 35 K, and is independent of the 
projection of selenium nuclear spin (mi = 0,±1/2). A 
comparison with P donors in silicon shows that Tj is be¬ 
tween one to two orders of magnitude longer for Se + in 
this temperature range. 

The temperature dependence of the electron spin T) 
has been well studied for shallow donors |T8], with differ- 


TLo = ui e S z - uj n I z + AS ■ I 
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Donor 

Raman temp¬ 
erature range (K) 

C (s^K 9 ) 
CT 9 

A (meV) Reference 

P 

2.6 - 6 

10 x 10 -8 

10.5 

m 

As 

5 - 11 

2 x 10“ 8 

19.7 

m 

Bi 

5 - 26 

6.6 x 10 -8 

33.8 

m 

Se+ 

5 - 35 

1.2 x 10 -8 

429 

this work 


TABLE I: Comparison of the strengths of the T 9 Raman elec¬ 
tron spin relaxation process for various shallow donors and 
Se+ in silicon, including the temperature range at which this 
mechanisms dominates Ti at X-band. Only a weak depen¬ 
dence is observed between the strength of the process, C, and 
the energy separation (A) between the ground ls(Ai) and 
first excited state (ls(T 2 )) of the donor. 


ent mechanisms dominating in different regimes of tem¬ 
perature and magnetic field. The Orbach two-phonon re¬ 
laxation process muni depends exponentially on the en¬ 
ergy gap, A = ls(T 2 ) — ls(Ai), between the ground state 
and first excited state of the donor: Ti a exp(A/fc s T). 
Therefore, though the Orbach process is the dominant T\ 
mechanism for phosphorus donors at X-band and at tem¬ 
peratures above ~8 K, the large value of A = 429 meV 
for Se + makes the Orbach process irrelevant here. A T 9 
Raman process dominates the phosphorus donor electron 
spin relaxation in the range 2 — 6 K at X-band j^L, aris¬ 
ing from two-phonon scattering via a continuum of ex¬ 
cited states. It has been observed in the shallow donors 
(see Table [T]) that this process has only a weak depen¬ 
dence on A, which dictates the spin mixing through a 
spin-orbit coupling. Our results on Se + are consistent 
with this observation - despite the much larger value 
of A, the strength of the T 9 Raman process is only 
marginally weaker than for the shallow donors. Neverthe¬ 
less, the electron spin T) for Se + in silicon remains longer 
than any of the shallow donors across the full tempera¬ 
ture range studied here, suggesting that it may also offer 
the longest spin coherence times. 

Electron spin coherence times (T 2 ) of ,7 Se + were mea¬ 
sured using a Hahn echo sequence (7t/2 — t — 7r — t — echo) 
and are shown in Figure [3] We find that T 2 is limited 
by spin relaxation {T\) for temperatures above ~12 K. 
Below this temperature, T 2 in these samples is limited 
by dipole coupling between Se + electron spins which are 
not refocussed by the Hahn echo sequence (a mechanism 
termed instantaneous diffusion (ID) (T2, [2T|). This effect 
can be mitigated by reducing the angle of the 8 2 pulse in 
the Hahn sequence (7r/2— t— d 2 —r—echo) ]12[lE2], provid¬ 
ing both a measure of the concentration of Se + electron 
spins, as well as a value of T 2 in the absence of ID. 
The concentration of Se + we obtain is 4 x 10 13 cm -3 , 
in good agreement with level of boron in this sample 
(noting that boron is required to ionise the selenium). 
The extended value of T 2 in the absence of ID, can be 
well described by a combination of spectral diffusion aris¬ 
ing from Tf-induced spin flips of neighbouring Se + spins, 



FIG. 3: Measured and extrapolated T 2 in 28 Si:' 7 Se + at 
B ~ 0.34 T. a) 1/T 2 as function of the refocusing pulse 
rotation angle # 2 . The intercept of the linear fit gives the 
extrapolated T 2 corresponding to the suppression of the in¬ 
stantaneous diffusion (ID), while the slope gives a Se + con¬ 
centration of 4 x 10 13 cm -3 , b) The value of T 2 is limited by 
Ti above ~12 K. Below this temperature, the extrapolated 
value of T 2 (in the absence of ID) can be well described by a 
combination of spectral diffusion arising from Ti-induced spin 
flips of neighbouring Se + spins, combined with a temperature- 
independent mechanism that limits T 2 to about 80 ms. 


combined with a temperature-independent mechanism 
that limits T 2 to about 80 ms. At the spin concentration 
used here, indirect flip-flop of neighbour electron spins 
is expected to limit T 2 to about ~1 s [12] . We there¬ 
fore postulate the temperature independent-mechanism 
observed arises from charge tunnelling in the sample, 
due to the high impurity concentration and compensa¬ 
tion in the material. Therefore, for lower donor con¬ 
centrations and using alternative means of ionising the 
Se, we anticipate the coherence times will be at least 
as long as those obtained for P (and indeed, longer, for 
the same temperature). The red solid line in Fig¬ 
ure [3] gives the temperature dependence of the extended 
value of T 2 (in the absence of ID). The model is ob¬ 
tained by taking into account the combined effects of 
the three mechanisms described above: i) spin relaxation 
of the central spin, yielding an echo decay of the form 
exp[— (2t/Ti)]; ii) a temperature-independent process of 
the form exp[— (2r/T 2i ii m )] with T 2i u m = 80 ms; and iii) 
a spectral diffusion process caused by spin relaxation of 
nearest neighbours, of the form exp[— (2t/Tsd) 2 ) with 
Tg D = KT 1 /[Se + ] and K = 7.6 x 10 12 s/cm 3 [Hl|21|. 


V. DONOR-ACCEPTOR RECOMBINATION 
DYNAMICS 

Pulsed ESR studies of donors at low temperatures of¬ 
ten employ optical excitation to generate free carriers, 










4 


# e- b) 

c) 


d) e - 

—5e + 

* 

—Se + 

4 pmf 

se°4f 


- fB° 

-if- 0 - 

B°—^— * 


O h + 
e) 




Af (s) 


FIG. 4: Schematic representation of the recombination mech¬ 
anisms in selenium doped silicon after above band gap illu¬ 
mination: (a) At the thermal equilibrium one of the electron 
from a selenium impurity moves to the acceptor B forming the 
DA pair (Se + B~). Under 1047 nm illumination electron-hole 
pairs are formed, (b) The excess of conduction carriers recom¬ 
bines through the impurities levels forming a non-equlibium 
state, (c) The process of DA recombination re-establish the 
thermal equilibrium, (d) Under 4 /im illumination an elec¬ 
tron from Se° is continuously pumped into the conduction 
band enhancing the rate at which the thermal equilibrium is 
re-established. The normalised concentration of selenium in 
the Sc + state ((Se + )) is monitored by plotting the normalised 
electron spin echo intensity, as a function of wait time, after 
a pulse of above band gap illumination, (e) in the dark, or (f) 
in the presence of 4 /im illumination. Red circles show nor¬ 
malised electron spin echo intensity and black lines are fits 
(see main text). 


as a way to ‘reset’ the electron spin to its thermal equi¬ 
librium state, allowing experiments to be repeated on 
timescales much faster than the intrinsic T\. However, 
for a compensated Se-doped sample, the free electron- 
hole pairs generated by the optical excitation are readily 
captured by the Se + and B“, and the ESR signal will only 
be recovered after charge equilibrium is reached through 


DA recombination. 

To study the dynamics of DA pairs, we measure 
the Se + electron spin echo as it recovers back to 
thermal equilibrium following above band-gap illumi¬ 
nation (1047 nm). This method has been previously 
used to study DA recombination of highly compensated 
phosphorus-doped silicon m showing that the recombi¬ 
nation time can be as long as ~ 10 3 s. Figure [ 4 ] gives a 
schematic representation of the DA recombination pro¬ 
cess after above band gap illumination. In this case, the 
system is brought out of equilibrium by the creation of 
free carriers via illumination (a). Due to the fast capture 
of such free carriers into the impurities levels, a non¬ 
equilibrium state is created after illumination (b). Ther¬ 
mal equilibrium is then re-established by DA recombina¬ 
tion (c). 

We first measured the photoconductivity directly fol¬ 
lowing illumination to obtain a photo-induced carrier life¬ 
time of 40 ms at 15 K (process (b) in Figure |4]). DA re¬ 
combination, a much slower process, was then measured 
by observing the electron spin echo intensity following the 
laser pulse, as shown in Figure |4j(e). A pulse of 100 mW 
for about 1 s was sufficient to ensure almost all the sele¬ 
nium and boron impurities were in the neutral state (i.e. 
zero initial spin echo signal). After ^1000 seconds fol¬ 
lowing the laser pulse, 60-70% of the echo intensity had 
recovered indicating that the majority of selenium-boron 
pairs had recombined. However, obtaining the full echo 
intensity required a waiting time of ~ 10 4 s — such non¬ 
exponential behaviour is consistent with a random dis¬ 
tribution of donor-acceptor nearest-neighbour distances. 

The problem of DA recombination has been solved an¬ 
alytically for shallow donors by Thomas et al. [23! and 
has been adapted here to selenium in the case of low com¬ 
pensation ([Se] 7g> [B]). In this case, the relative distance 
between donors and acceptors is much larger than the 
Bohr radius of both electron and hole bound the to impu¬ 
rities and the wavefunctions in the ground state may be 
considered unperturbed. For instance, for concentrations 
of ~ 10 15 cm -3 and low compensation, the most prob¬ 
able distance between DA nearest-neighbours is about 
50 nm which is much larger than a* H = 4.2 nm, the ef¬ 
fective Bohr radius of the bound hole in silicon. Under 
this approximation, the calculation of the recombination 
rate ( W ) for each individual DA pair only involves the 
optical matrix element M between the wavefunctions of 
the donor and acceptor in the neutral condition and we 
have W(r) = \M\ 2 . The optical matrix element can be 
then taken to be proportional to the value of the hole 
wavefunction at the donor site, giving 

W(r) = Wo exp , (2) 

l a H J 

where Wq has been determined to be between 10 5 s _1 
and 10 3 s" - 1 m- The effect of an ensemble of selenium 
impurities surrounding the acceptor can be seen as set of 
independent recombination channels with a certain dis¬ 
tribution of recombination rates. By integrating over 
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a random distribution of donors, the dynamics can be 
shown to follow: 


<B°(t)) 

[B] 


= exp 


4^[S e] \e~ w ^ r)t 


*} 



(3) 


(Se+(t)) _ (B°(t)) 

[B] [B] 


(4) 


where (Se + (f)) and (B°(f)) are the concentrations of Se + 
and B _ respectively. We simulated the dynamics of the 
recombination process using the equations above, tak¬ 
ing the known concentration of selenium impurities, and 
find a good fit to the experimental data. The value of 
the relaxation rate Wq used in the fit was 10 4 s _1 , at 
15 K, which is consistent with previous values observed 
for similar samples doped with shallow donors m 
The approach to thermal equilibrium (all boron ionised 
and an equal number of Se + ) can be accelerated over that 
resulting from DA pair recombination by using infrared 
radiation at about 4 /j m: an electron is continuously ex¬ 
cited from the neutral selenium into the conduction band, 
whereupon it either returns to re-neutralise the selenium, 
or is captured by an acceptor (Figure |4](d)). We find 
that under weak (10 /tW) illumination at 4 /mi, the rate 
at which equilibrium charge conditions are established is 
two orders of magnitude faster then the direct DA re¬ 
combination. Within the first second after the 1047 nm 
flash, about 50 % of the echo signal is recovered and 
thermal equilibrium state is reached after about 100 s. 
The recombination is well fit by a bi-exponential decay 
with Ti = 1.6 s and r 2 = 40 s. The effect of 4 /tm il¬ 
lumination is further evidence that the slow observed 
recombination dynamics are limited by long-lived (i.e. 
distant) DA pairs. Moreover the illumination sequence 
(1047 nm + 4 /tm) could be used in order to study the 
nuclear spin properties of neutral selenium, ,7 Se°, fol¬ 
lowing analogous methods to those used to study ionised 
P nuclear spins through electrically detected magnetic 
resonance [53]. In shallow donors such as P, the hyper- 
fine interaction with the electron spin strongly limits the 
nuclear spin coherence, such that nuclear spin coherence 
times for P + are strongly enhanced [7j. Similarly, one 
could expect the nuclear spin coherence time for neutral 
selenium (in its singlet ground state) to be comparably 
long. 


VI. CONCLUSIONS 

We have found that the electron spin relaxation times 
of Se + are the longest observed for donors in 28-silicon 


- this is explicitly demonstrated for temperatures in 
the range 5-30 K, and is expected to hold true for 
lower temperatures. The energy splitting, A, between 
the ground and first-excited valley state is an order 
of magnitude larger in Se + than for Group V donors 
in silicon — this effectively removes the Orbach spin 
relaxation mechanism, though its effect in reducing 
the strength of the Raman T 9 process is rather weak, 
providing motivation for further theoretical work to 
understand this quantitatively. As a result of this 
increase in T), the electron spin coherence times in Se + 
are shown to be significantly longer than for Group V 
donors for temperatures above 10 K, and this trend is 
likely to extend to lower temperatures for samples with 
lower Se concentrations than those studied here. The 
large value of A should also result in a dramatically 
reduced Stark shift [25] compared to shallow donors, 
such that the electron and nuclear spin coherence of 
Se + could be relatively unaffected by charge noise in 
nano-devices, despite the large value of the hyperfine 
coupling. It is also possible to spectrally resolve the 
ground state hyperfine coupling in 28 Si: 77 Se + using the 
1s(T 2 )T 7 absorption transition [2], so if a tuneable 
source were available at ~ 2.9 /tm this could be used for 
fast and efficient hyperpolarization of both the electron 
and nuclear spin. Finally, we have investigated optical 
methods to manipulate the donor charge state (from 
Se + to Se°), examining this through donor-acceptor 
recombination following above band gap illumination, 
and showing how Se° can be ionised using 4 /tm illu¬ 
mination. This suggests a possible mechanism to first 
bring selenium donors into an ESR-active Se + state, 
and subsequently neutralise them leaving a potentially 
long-lived ,7 Se nuclear spin. 
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